Electric field scaling at a B = metal-insulator transition in two dimensions 
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The non-linear (electric field-dependent) resistivity of the 2D electron system in silicon exhibits 
scaling as a function of electric held and electron density in both the metallic and insulating phases, 
providing further evidence for a true metal-insulator transition in this 2D system at B = 0. Compar- 
ison with the temperature scaling yields separate determinations of the correlation length exponent, 
v ~ 1.5, and the dynamical exponent, z ~ 0.8, close to the theoretical value z = 1. 
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Conventional wisdom holds that a two-dimensional 
system of noninteracting electrons is always localized in 
zero magnetic field at zero temperature no matter how 
weak the disorder jl). If so, there should be no metal- 
insulator phase transition at B — in an infinite sam- 
ple in 2D. However, recent experiments m have shown 
clear signatures of a metal-insulator transition in a high- 
quality two-dimensional electron system (2DES) in sili- 
con in zero magnetic field. The resistivity, p, was found 
to scale with temperature, so that 

p(T,n s ) = ft (\5 n \/T b ) with b = 1/zv. (1) 

Here n s is the electron density, 5 n = (n s — n c )/n c , n c is 
the critical electron density at the metal-insulator tran- 
sition, z is the dynamical exponent, and v is the cor- 
relation length exponent. On the insulating side of the 
transition, the resistivity was found to have a temper- 
ature dependence characteristic of hopping in the pres- 
ence of a Coulomb gap in the density of states j| due to 
strong electron-electron interactions. Such scaling behav- 
ior is characteristic of a true phase transition and implies 
that there is a diverging correlation length, £ ~ |(5 n | _ly , 
at the transition. The behavior strongly resembles the 
superconductor-insulator phase transition in thin disor- 
dered films m, as wen as the phase transition between 
quantum Hall liquid and insulator [j?]Je]j9|JlO| . These re- 
sults are in apparent contradiction with the scaling the- 
ory [111 which, however, ignores electron-electron inter- 
actions. In the 2DES in silicon, electron-electron inter- 
actions provide the dominant energy, E ee , at low elec- 
tron densities: at n s ~ 10 11 cm~ 2 , E ee ~ 5 meV while 
the Fermi energy is only ~ 0.6 meV. For comparison, 
in GaAs/AlGaAs heterostructures these energies are ap- 
proximately equal to each other at the same electron den- 
sity: E ee w Ep w 3.5 meV. We note that recent theo- 
retical studies [ pd]p^ ] have shown that strong electron- 
electron interactions can cause derealization. 

The resistivities reported in Refs. ||||] were obtained 
in the linear regime, i.e., in the limit of zero electric field, 



E —y 0. When the electric field is strong, however, the 
effective temperature of the electrons becomes different 
from the lattice temperature. A general scaling analysis 
of the Coulomb interacting quantum critical point |l5| ] 
shows that the resistivity should also scale with the elec- 
tric field, but the scaling exponent in this case will be 
l/[(z+ (instead of 1/zv in the case of temperature 
scaling) so that 

p(E,n s ) = f 2 {\5 n \/E a ) with o = l/[(«+ (2) 

Knowing both zv and (z + l)i/, one should be able to de- 
termine the exponents z and v separately. It is worth not- 
ing that this was done recently for the superconductor- 
insulator transition in thin disordered films Q. 

In this Letter we report measurements of the nonlin- 
ear resistivity of a two-dimensional electron system as 
a function of electric field at B — 0. Scaling is ob- 
served as a function of electric field and electron den- 
sity about a critical point which corresponds to the same 
critical density about which temperature scaling was re- 
ported earlier By comparing the electric field scal- 
ing with the temperature scaling, separate determina- 
tions (rather than their product) were obtained for the 
correlation length exponent, v « 1.5, and the dynamical 
exponent, z » 0.8, close to the theoretical value z = 1 
for strongly interacting Coulomb systems. At the tran- 
sition (n s = n c ), the resistivity is independent of electric 
field and close to 3/i/e 2 . These data provide additional 
strong evidence for a true metal-insulator transition in 
high-mobility 2DES in silicon in the absence of a mag- 
netic field. 

The samples used in this work were high-mobility sil- 
icon metal-oxide-semiconductor field-effect transistors 
(MOSFET's) with maximum electron mobility p max rj 
35 - 40, 000 cm 2 /Vs similar to those used in Refs. ||,§. 
The resistivity was determined as a function of electric 
field by measuring current versus voltage (7 — V curves) 
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FIG. 1. Nonlinear current- voltage curves at several electron densi- 
ties around the metal-insulator transition. T = 0.22 K. 

for many values of n s using a standard four-terminal dc 
technique. Typical curves are shown in Figure |. The 
middle, solid curve is close to linear and separates all 
I — V curves into two groups exhibiting different types 
of nonlinear behavior. Note that there is a remarkable 
symmetry about the middle curve. This symmetry is 
reminiscent of that observed for the I — V curves in the 
vicinity of the QHE-insulator phase transition and 
is consistent with the existence of a critical point at a 
critical electron density n c . In Ref. P], the symmetry 
was attributed to charge-flux duality. In our case, its 
physical origin is unclear and merits theoretical study. 

Curves of resistivity p vs electric field E for different 
electron densities are shown in Fig. [2] (not all curves 
measured are shown in order to avoid too high a den- 
sity of points). Here, the resistivity is determined from 
p = (V/I) ■ (W/L) (W is the sample width and L is 
the distance between potential contacts) and the electric 
field E = V/L. Again, there is clearly a critical elec- 
tron density which separates two distinct density regions 
characterized by different types of resistivity behavior 
as a function of electric field. At the critical point, the 
resistivity (solid curve in Fig. [j], closed circles in Fig. ||) 
is virtually independent of the electric field and close 
to 3/i/e 2 . All curves below this line are characterized 
by dpj 'dE > while all curves above are character- 
ized by dp/dE < 0. The observed critical behavior of 
p as a function of E is very similar to the critical be- 
havior of the resistivity with temperature observed near 
the metal-insulator transition in 2DES in silicon [^|J|] 
(with the only difference that the resistivity at the crit- 
ical point is _E-independent but not T-independent at 
T > 2 K). It also resembles the behavior found near the 
superconductor-insulator phase transition in thin metal- 
lic films H and the quantum Hall effect-insulator (QHE- 
I) transition in GaAs/AlGaAs heterostructures [MMflOj . 
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FIG. 2. Resistivity as a function of electric field at Tl s = 7.81, 

7.92, 8.03, 8.14, 8.25, 8.36, 8.47, 8.70, 8.91, 9.13, 9.35, 9.57, 
9.79, 10.34, and 10.78 x 10 10 m r 2 at T = 0.22 K. 

A remarkable property is that plotting the resistiv- 
ity against the scaling variable, \5 n \/E a , causes all the 
curves to collapse onto two distinct branches, as shown 
in Fig. H (a). Note that points for E — ► 0, where p sat- 
urates due to finite temperature, have been omitted. As 
expected, the scaling fails at electron densities far from 
the critical point where the system is no longer in the 
critical regime, as seen on the lower curve in Fig. || (a). 
The upper branch corresponds to n s < n c and the lower 
one to n s > n c . At the transition (n s — n c ), the re- 
sistivity is close to 3h/e 2 . It is interesting to note that 
almost the same value was recently reported in Ref. |l4| ] 
for a transition between weak and strong localization in 
disordered metallic films. The exponent a was varied to 
obtain the best visual collapse, yielding a — 0.37 ± 0.01. 

Repeating the procedure reported in Refs. H[|, the 
temperature dependence of the resistivity measured in 
the linear regime (E — > 0) was used to obtain scaled 
curves of the resistivity as a function of the temperature- 
dependent scaling variable, \5 n \/T b . This is shown in 
Fig. U (b), where it is evident that all curves again col- 
lapse onto two distinct branches, confirming the earlier 
results. We note that the data scale well only at temper- 
atures below 1 K, presumably because at higher temper- 
atures the system is outside the critical regime. The best 
collapse was achieved for b = 0.83 ± 0.08. 

As mentioned earlier, these two scaling analyses allow 
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one to determine both the dynamical exponent and the 
correlation length exponent. The values obtained in 
our experiments are v = 1.5 ± 0.1 and z = 0.8 ± 0.1. 
While the correlation length exponent has not been the- 
oretically predicted for this system, the value of the 
dynamical exponent in a strongly interacting 2D system 
was predicted to be z — 1 |]l3), close to what we ob- 
serve. In the majority of interacting 2D systems where 
it has been determined, the exponent z has also been 
found to be close to 1 (see, e.g., Ref. Q for the case of a 
superconductor-insulator transition and Ref. |l5| for the 
transition between two neighboring QHE plateaux). The 
same result for the case of hopping conductance follows 
from theory of Polyakov and Shklovskii Jll| who found 
the ratio of scaling exponents a/b — 1/2 (see Eqs. [j] and 
||) which is equivalent to z = 1. 

We note that the above considerations were based on 
the assumption that there is a diverging correlation time 
at the critical point which controls the dynamics Ej) . 
An alternative picture involving hydrodynamic heating 
was developed in Ref. |lj} which assumes that the elec- 
tron temperature is governed by phonon emission. In the 
case of silicon, this model gives a/b = 1/3 rather than 
a/b = 1/2. Some influence of phonons may be responsi- 
ble for the ratio a/b somewhat less than 1/2 (a/b « 0.45) 
found in the present work. A crossover between the two 
regimes will be reported elsewhere jl8| . 

In summary, we have shown that in the absence of a 
magnetic field, the resistivity of the 2DES in high-quality 
silicon MOSFETS scales as a function of electric field 
and electron density, exhibiting critical behavior about 
the same point as the temperature scaling. At the crit- 
ical point the resistivity is independent of electric field 
and close to 3h/e 2 . Comparison of scaling of the resistiv- 
ity with temperature and scaling with electric field yields 
separate determinations of the correlation length expo- 
nent, v « 1.5, and the dynamical exponent, z w 0.8, 
close to the theoretical value z = 1. These results pro- 
vide additional strong evidence of a true metal-insulator 
phase transition in a high-quality 2DES in silicon in zero 
magnetic field. 
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handling and preparation of the manuscript. This work 
was supported by the US Department of Energy under 
Grant No. DE-FG02-84ER45153. 



FIG. 3. (a) Demonstrating scaling with electric field, the resis- 
tivity at 0.22 K is plotted as a function of \5 n \/E a for a = 
l/[(z + 1)^"] = 0.37. (b) Demonstrating scaling with temper- 
ature, the linear resistivity (E — > 0) is shown as a function of 
\8 n \/T^ for b = \/zv = 0.83. Electron densities are in the 
range 7.81 to 10.78 X 10 10 cm~ 2 . 



[1] E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. 

V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 
[2] S. V. Kravchenko, G. V. Kravchenko, J. E. Furneaux, V. 



3 



M. Pudalov, and M. D'lorio, Phys. Rev. B 50, 8039 (1994). 

[3] S. V. Kravchonko, W. E. Mason, G. E. Bowker, J. E. 
Furneaux, V. M. Pudalov, and M. D'lorio, Phys. Rev. B 
51, 7038 (1995). 

[4] A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 (1975). 

[5] Y. Liu, K. A. McGreer, B. Nease, D. B. Haviland, G. Mar- 
tinez, J. W. Halley, and A. M. Goldman, Phys. Rev. Lett. 
67, 2068 (1991). 

[6] A. Yazdani and A. Kapitulnik, Phys. Rev. Lett. 74, 3037 
(1995). 

[7] T. Wang, K. P. Clark, G. F. Spencer, A. M. Mack, and W. 

P. Kirk, Phys. Rev. Lett. 72, 709 (1994). 
[8] D. Shahar, D. C. Tsui, M. Shayegan, R. N. Bhatt, and J. 

E. Cunningham, Phys. Rev. Lett. 74, 4511 (1995). 
[9] D. Shahar, D. C. Tsui, M. Shayegan, E. Shimshoni, and S. 

L. Sondhi (unpublished). 
[10] L. W. Wong, H. W. Jiang, N. Trivedi, and E. Palm, Phys. 



Rev. B 51, 18 033 (1995). 
[11] D. L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (1994). 
[12] A. L. Efros and F. G. Pikus, Solid State Commun. 96, 183 

(1995). 

[13] S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, to 
be published. 

[14] S.-Y. Hsu and J. M. Valles, Jr., Phys. Rev. Lett. 74, 2331 
(1995). 

[15] H. P. Wei, L. W. Engel, and D. C. Tsui, Phys. Rev. B 50, 
14 609 (1994). 

[16] D. G. Polyakov and B. I. Shklovskii, Phys. Rev. B 48, 
11 167 (1993). 

[17] S. M. Girvin, private communication; see also E. Chow, H. 
P. Wei, S. M. Girvin, and M. Shayegan, Phys. Rev. Lett., 
in press. 

[18] D. Simonian, S. V. Kravchenko, and M. P. Sarachik, in 
preparation. 



4 



